Search for Neutral Heavy Leptons in a High-Energy Neutrino Beam 
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A search for neutral heavy leptons (NHLs) has been performed using an instrumented decay 
channel at the NuTeV (E-815) experiment at Fermilab. The decay channel was composed of helium 
bags interspersed with drift chambers, and was used in conjunction with the NuTeV neutrino detec- 
tor to search for NHL decays. The data were examined for NHLs decaying into muonic final states 
{^J.^J.l', fiev, fin, and up); no evidence has been found for NHLs in the 0.25 - 2.0 GeV mass range. 
This analysis places limits on the mixing of NHLs with standard light neutrinos at a level up to an 
order of magnitude more restrictive than previous search limits in this mass range. 

PACS numbers: 13.15.+g,13.35.Hb,14.60.Pq 
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Various extensions to the Standard Model predict 
neutral heavy leptons (NHLs) which can mix with the 
standard light neutrinos. In these extensions, the NHLs 
are weak isosinglets that do not couple directly to the 
Z and W bosons, but can decay via mixing with the 
Standard Model neutrinos. Figure |l| shows one possible 
set of tree-level diagrams for NHL production and decay. 
The NuTeV (E815) neutrino experiment at Fermilab has 
made a sensitive search for these NHLs by combining the 
capabilities of a high intensity neutrino source with an 
instrumented decay region. 

In these extended models [0 , the NHL lifetime depends 
on the mixing parameter |J7p and the mass of the NHL. 
They are expected to decay (e.g. Fig. [|b) into a neutrino 
and two charged leptons, into a lepton and two quarks, 
or into three neutrinos. 

NHLs may be created in the NuTeV beamline by the 
decays of secondary mesons produced by the Tevatron 
proton beam. During the 1996-1997 fixed-target run at 
Fermilab, NuTeV received 2.54 x 10^^ 800 GeV protons 
on a beryllium oxide production target with the detector 
configured for this search. A sign-selected quadrupole 
train focused secondary tt and K mesons down a beam- 
line 7.8 mrad from the primary proton beam direction. 
1.13 X 10^* protons were received with the magnets set 
to focus positive mesons, and 1.41 x 10^* protons with 
negative meson focusing. The mesons could decay into 
NHLs as shown in Fig. The production of secondary 
pions and kaons was simulated using the parameteriza- 
tion in 1^; the Decay Turtle program [Q simulated the 
propagation of charged particles through the beamline. 

NHLs may also be produced by prompt decays of 
charmed mesons produced by incident protons on the 
BeO target and proton dumps. These processes were 
simulated using a Monte Carlo program based on mea- 
sured production cross sections [g| . The effects of decay 
phase space, NHL polarization, and helicity suppression 
[|6[ were included in the simulation of the production and 
decays. For NHLs of mass 1.45 GeV from D meson decay, 
the average momentum was ^100 GeV; for the 0.35 GeV 
NHLs coming mainly from K decay, the average momen- 
tum was ~140 GeV. 

This analysis reports the results of a search for NHLs 
with masses between 0.25 to 2.0 GeV, which decay with a 
muon in the final state. The primary NHL decay modes 
of this type are /iei^, /i/ii', /ivr, and [ip. In the standard 
mixing model 0J^, the expected ratio of decay rates for 
/iTT : [lev : [ipv is 3.5 : 1.6 : 1.0. The \ip channel is 
significant only for NHL masses above 1 GeV, where it 
is comparable to the /i7r mode. 

In order to detect NHL decays, an instrumented de- 
cay region (the "decay channel") was constructed 1.4 km 
downstream of the production target. The decay channel 
(Fig. ^ was designed to identify signatures characteristic 
of NHLs: a neutral particle decaying to two charged par- 
ticles within the decay channel. A 4.6 m x 4.6 m array 



of plastic scintillators vetoed charged particles entering 
from upstream. The decay channel was 34 m long, in- 
terspersed with 3 m X 3 m multi-wire argon-ethane drift 
chambers positioned at 15.4 m, 25.2 m, 35.5 m, 36.8 m, 
and 38.5 m downstream of the veto array. Tracks were 
reconstructed from drift chamber hits and grouped to 
form candidate decay vertices. Typically, a vertex from 
a NHL of mass 1.15 GeV was reconstructed with a reso- 
lution of 0.04 cm in the transverse direction and 28 cm in 
the longitudinal direction. The region between the drift 
chambers was filled with helium contained in cylindrical 
plastic bags 4.6 m in diameter. By displacing the air with 
helium, the number of background neutrino interactions 
between the chambers was reduced by a factor of 7. 

The Lab E neutrino detector located immediately 
downstream of the decay channel, provided final state 
particle identification, energy measurement, and trigger- 
ing. This detector consisted of a 690-ton iron-scintillator 
sampling calorimeter followed by a toroidal muon spec- 
trometer. Drift chambers were positioned every 20 cm 
along the length of the calorimeter, and 2.5 cm-thick liq- 
uid scintillator counters were interleaved with 10 cm steel 
plates along the length of the calorimeter. The trigger se- 
lected events with a penetrating muon or at least 6 GeV 
of hadronic or electromagnetic energy in the calorimeter. 
Offline, hits in the calorimeter drift chambers were ana- 
lyzed to determine if they formed a track (consistent with 
a muon) or a cluster of hits (consistent with an electron 
or pion shower). 

Tracks found in the decay channel were linked to clus- 
ters or tracks in the calorimeter. Pulse height infor- 
mation from the scintillation counters was used to de- 
termine hadronic or electromagnetic energy deposition; 
the hadronic energy resolution of the calorimeter was 

= 0.024 and the electromagnetic energy res- 
olution was = 0.040 Muon energies were de- 
termined either by the toroid spectrometer (with 11% 
resolution) or by range (with a resolution of 155 MeV). 
The energy of muons which exited through the side of 
the calorimeter was determined from the track's multi- 
ple scattering, with a resolution of 25% for energies lower 
than 85 GeV. 

The following processes were the main backgrounds 
to this search: 1) neutrino interactions within the he- 
lium; 2) neutrino interactions in the material upstream 
of the decay channel, where a neutral kaon survived to 
the decay channel and decayed in the fiducial volume; 
and 3) neutrino interactions in the material surround- 
ing the decay channel. Neutrino interaction backgrounds 
could arise from low multiplicity resonance production 
or from deep inelastic scattering (DIS). Resonance pro- 
duction was simulated using the calculations of Belusevic 
and Rein [^; such events were characterized by a high- 
energy muon accompanied by a low-energy pion track. 
The Lund Monte Carlo program was used to simulate 
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DIS events [|T^ . While these were generally of high mul- 
tiplicity, events with few or poorly-reconstructed tracks 
could contribute to the two-track background. The de- 
tector was simulated using the GEANT Monte Carlo, 
which produced a hit-level simulation of raw data includ- 
ing beam and cosmic ray related noise hits. These were 
processed using the same analysis routines as those used 
for the data. 

Backgrounds from cosmic rays were estimated using 
a sample of cosmic ray muons which interacted in the 
target calorimeter; an upper limit of 10~^ background 
events was determined. We have therefore ignored cosmic 
rays in the final background estimate. 

Event selection criteria maintained high efficiency for 
the NHL signal while minimizing known backgrounds. 
Cuts fell into two broad categories: reconstruction and 
kinematic. 

Reconstruction cuts isolated events with a two-track 
decay vertex within the decay channel fiducial volume 
and no charged particle identified in the upstream veto 
system. The two tracks were required to be well- 
reconstructed, have an accompanying energy measure- 
ment, and form a common vertex. Large angle tracks 
arising from cosmic rays were removed by requiring 
tracks to form an angle of less than 0.1 radians with 
the beam direction. Exactly two tracks were required, 
both projecting to the calorimeter, with at least one of 
the two identified as a muon. To ensure good particle 
identification and energy measurement, muons were re- 
quired to have an energy greater than 2.0 GeV; an en- 
ergy greater than 10.0 GeV was required for electrons 
or hadrons. (The latter cut eliminated the low-energy 
pions associated with resonance production.) The recon- 
structed two-track decay vertex was required to be at 
least 1 m from any material in the drift chambers. 

The kinematic cuts were designed to remove the re- 
maining DIS and resonance backgrounds. The effective 
scaling variables Xcs and WcS were calculated for each 
event under the following assumptions: 1) the event was 
a neutrino charged current interaction {vN ^N'X); 
2) the highest energy identified muon was the outgoing 
particle from the lepton vertex; and 3) the missing trans- 
verse momentum in the event was carried by an unde- 

tected final state nucleon. Specifically, XcS = 2m "i^" 

and Wcs = + 2TOpi^vis - Qvis' where Ql^^ is the re- 

constructed momentum transfer squared, uip is the mass 
of the proton, and i^vis is the reconstructed hadron en- 
ergy. Requiring a;cff < 0.1 and reduced backgrounds 
from DIS; requiring WoS > 2.0 GeV removed quasielestic 
and resonance backgrounds. Since we could not recon- 
struct the true mass of the NHL due to the missing 
neutrino, a cut was applied on the "transverse mass," 
Tot = \pt\ + v^p^ + TOy. pt is the component of the 
total momentum perpendicular to the beam direction, 
and m-v is the invariant mass for the two charged tracks. 



Requiring tot < 3.0 GeV removed additional DIS back- 
ground. 

The NHL and background acceptances and reconstruc- 
tion efficiencies were calculated from the hit-level Monte 
Carlo simulation. For a range of NHL masses from 0.25 
to 2.0 GeV, mesons were allowed to decay in the beam- 
line to muons and NHLs. The NHLs from these decays 
were propagated to the decay channel, weighted by their 
decay and the meson production probability, and allowed 
to decay in the decay channel producing hits in the drift 
chambers and energy deposition in the Lab E calorime- 
ter scintillation counters. These simulated raw data were 
reconstructed with the same programs used for the data 
including the above reconstruction and kinematic cuts. 
For NHL decays generated within the decay channel fidu- 
cial volume, the average acceptance was 23% for masses 
between 0.25 and 2.0 GeV. From the background simula- 
tion, the calculated background contributions for the full 
data sample after cuts are listed in Table |[ 

Neutrino events in the drift chambers and DIS (>2 
track) events in the helium have been used to cross check 
the Monte Carlo calculations of backgrounds and efh- 
ciencies. For example, as a check on the two-track recon- 
stuction efficiency, we have compared the fourteen two- 
track data events with vertices in the drift chambers to 
the Monte Carlo prediction of 15 ± 2 events. Figure ^ 
shows comparisons of distributions for the previously in- 
troduced DIS variables, a^cff, W^off , Qls as well as tot and 
ycs = ^ • Data events with > 2 tracks also showed 
good agreement with the Monte Carlo prediction: 280 
events were predicted, and 275 were observed. 

For the full NuTeV data sample, no data events 
passed all cuts. This is consistent with the expected 
background of 0.57 ± 0.15 events. We have set limits 
from this null result by using the Monte Carlo predic- 
tion for Npy^edimj^Hh, |C^P)j the normalized number of 
NHL events expected in the decay channel as a func- 
tion of mass and |C/p. The 90% confidence level limit 
for a given mass was calculated for the null observa- 
tion by finding the |f/|^ value for each mass such that 
A^piod('7iNHL, = 2.3 events. 

The statistical 90% confidence level limit was modi- 
fied by the addition of systematic uncertainties on the 
number of expected NHLs. These uncertainties are sum- 
marized in Table |l^. Experimental uncertainties on the 
D production cross section and on the K meson fiux 
Monte Carlo Q| were the dominant contributions to the 
systematic error. The systematic uncertainties were in- 
corporated into the 90% confidence level limit by adding 
in quadrature a fractional error term corresponding to 
the systematic uncertainty for a given mass. Since the 
NHL rate was proportional to \U\'^, adding the system- 
atic uncertainties increased the \U\'^ limit by 4%(14%) at 
0.35(1.45) GeV. 

Figure Q shows the limits obtained from this search for 
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the NHL-j/^ mixing parameter jiJ^I as a function of the 
mass of the NHL. The hmits are for NHL decay modes 
containing a muon in the standard mixing model pre- 
sented above. Also shown in Fig. ^ are the results of 
previous experiments. |p^p^ 

In summary, we have made a sensitive search for NHLs 
in the mass range from 0.25 to 2.0 GeV that mix with 
muon neutrinos. No evidence for NHL production was 
observed. New limits have been set that are up to an 
order of magnitude better than previous searches in the 
lower mass range. 
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Energy and the National Science Foundation. We thank 
the staff of FNAL for their contributions to the construc- 
tion and support of this experiment during the 1996-97 
fixed target run. 
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background source description 


number of events 
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K*^ punch-through 
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0.56 ±0.15 
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0.002 ± 0.001 
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FIG. 1. Feynman diagrams for the (a) production and (b) 
decay of neutral heavy leptons (L^). |f/|^ is the mixing be- 
tween the muon neutrino and the NHL. (The alternative al- 
lowed decay via a Z boson is not shown.) 



NC 

Beam 

line 



Veto 
Wall 



Drift 
Chamber 



Drift 
Chamber 



Drift 
Chambers 



Helium Bag 



Helium Bag Lab E Deteetor 



FIG. 2. A schematic diagram of the NuTeV decay channel. 
The beam enters from the left; at the far right is the Lab 
E calorimeter. An example NHL decay to fj,Tr is also shown. 
The event appears as two tracks in the decay channel, and a 
long muon track and a hadronic shower in the calorimeter. 




GeV/rJ 



0' 



Z posit 



25 
meter 



0.02 0.04 0.06 0.08 0.' 



FIG. 3. Kinematic distributions for data (crosses) and MC 
DIS background (dashed) events with two-track vertices re- 
constructed in the decay channel drift chambers. MC events 
are absolutely normalized to the number of protons on target. 
The DIS reconstructed variables j/eff, WbS, Qeffi ^eff and mr 
are defined in the text. Z positions are referenced to the veto 
array; spikes in the Z distribution correspond to the locations 
of the decay channel drift chambers. 
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FIG. 4 

ing of NHLs to Standard Model left-handed muon neutrinos 
as a function of NHL mass. The solid line in the figure corre- 
sponds to the limit with systematic errors included. 
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